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The collective phase response to a macroscopic external perturbation of a population of interacting 

nonlinear elements exhibiting collective oscillations is formulated for the case of globally-coupled 

oscillators. The macroscopic phase sensitivity is derived from the microscopic phase sensitivity of 

QQ ' the constituent oscillators by a two-step phase reduction. We apply this result to quantify the 

C^ ' stability of the macroscopic common-noise induced synchronization of two uncoupled populations 

f^ , of oscillators undergoing coherent collective oscillations. 

PACS numbers: 05.45.Xt, 82.40.Bj, 05.40.Ca 

QQ ■ Nature possesses various examples of systems coinposed of many nonlinear elements in which, through mutual 

interactions, coherent collective dynamics emerges [l|, y, S \% la, m, LZl, Isl, Isl, llO, lul, llJ, ll3|. The collective dynamics of 



^^ the population is generally not a simple superposition of the microscopic dynamics of the constituent elements; due to 

v^ nonlinearities, it can differ qualitatively. Understanding how the microscopic properties give rise to the macroscopic 
collective properties of the system, and formulating the theory in a simple closed form at the macroscopic level is a 
fundamental problem in nonlinear dynamics. 

If the population undergoes a stable macroscopic limit-cycle oscillation, it can most simply be described by a single 
macroscopic phase variable. By knowing the collective phase response of the macroscopic oscillation with respect to 
external perturbations, we can describe the dynamics of the population by the phase picture, much like for single 
microscopic oscillators [l|, y, y, |j, |5| . The collective phase response can be measured operatively by macroscopically 
applying an external perturbation to the population without recourse to the microscopic details [l|, ll2|, ll3| . However, 
l/~j ' to understand how the microscopic dynamics of the constituent elements and their mutual interactions conspire to 
OO . become a macroscopic phase response, it is necessarily to develop a statistical mechanical approach. 
f^ ' In this paper, we formulate this problem for globally-coupled oscillators as a simple model of macroscopic popu- 

. . lations exhibiting collective oscillations. Global coupling is the simplest form of mutual interactions, which is both 
C ■ conceptually and practically important because it serves as a starting point for theoretical analysis of various network- 
^^ . coupled dynamical systems and because it is realized in experimental systems [l|, l3, ISl U, IS, ISl ■ We quantitatively 
derive the collective phase sensitivity, which is the linear response coefficient of the macroscopic collective phase re- 
sponse, from the microscopic phase sensitivity of the constituent oscillators using a two-step phase-reduction on the 
micro- and macroscopic scales. 

The collective phase sensitivity is a fundamental quantity characterizing the response of a macroscopic limit- 

^ ■ cycle oscillation to weak external perturbations. As one application, we analyze macroscopic common-noise induced 

synchronization of the collective oscillations between two uncoupled populations. Though not treated in this paper, 

the entrainment of a population to external forcing or synchronization between populations can also easily be analyzed 

using the collective phase sensitivity. 

Let us give a general definition of the collective phase sensitivity first. We consider a population of globally-coupled 
noisy identical limit-cycle oscillators collectively exhibiting stable macroscopic limit-cycle oscillations. We assume that 
the effect of coupling, noise, and perturbations is sufficiently weak, so that the orbit of each microscopic oscillator 
is always near its limit cycle. We can then describe each oscillator by the phase 0, and the whole population by a 
single-oscillator phase probability density function (PDF) P{4), t) of the constituent oscillators. As the population 
undergoes limit-cycle oscillations, P{(t), t) is given in the form of a rotating wave packet of constant shape (see Figs. [T] 
and [3]), 

p(0,i) = /(^-<i>(i)), <^{t) = nt + Q, (1) 

where f{(f) represents the wave packet, <^{t) the collective phase of the population (location of the wave packet) at 
time t, Vt the frequency of the collective oscillation, and 8 the initial collective phase. We assume e [— 7r,7r] and 
impose a periodic boundary condition f{4> + 27r) = f{(j)). 



At the instant that the collective phase is $0j we apply a macroscopic perturbation of magnitude and direction s 
that acts uniformly on all oscillators. The shape of P((/), t) transiently deforms, but due to the stability of the collective 
oscillations, eventually settles into its stationary shape, Eq. ^. We denote the perturbed PDF after a sufficiently 
long time ip as Pi{4>^t + to) = / ('/' ^ *&i(^Oj •Sj^o)); and the unperturbed PDF after time to as P2{4>,t + to) — 
f {4> — $2(*&07 ^o)), where $2 = $o + ^^o- The collective phase response function G'($o; •s) is defined as the asymptotic 
difference of $i($o, «, ^o) and $2(^0, io): G($o;s) = limjo^oo [$i(*o, «, ^o) - *2(*o,io)] [I- This is in general a 
function of $0 and s. When the impulse magnitude \s\ is small enough, the collective phase response function is 
proportional to s, and is given by G(<I>o; s) = C(^o) • S; where C(^o) is the collective phase sensitivity [2]. 

We consider globally-coupled oscillators described by the following Langevin equation: 



1 ^ , 

X,{t) = F{X,{t)) + - ^ V{Xu{t)) + V^V.it) + sit), 



(2) 



k=l 



for j = 1, • • • , iV, where Xj is the state of the jth oscillator, F{X) the dynamics of a single oscillator, V{X) the 
contribution from each oscillator to the mean field, Dg the noise intensity, rjj (t) an independent zero-mean Gaussian 
white noise of unit intensity, which satisfies {'nj{t)'nk{t')) = Mj,fc^(t — t') where I is the unit matrix, and s{t) the 
macroscopic external perturbation common to all oscillators. When isolated, each oscillator has a stable limit-cycle 
Xo{t + T) ^ Xo{t) with period T — 2tt/uj. We define a phase 4>{X) E [— tt, tt] around the limit cycle, which increases 
with a constant natural frequency uj l^. We assume that the effect of coupling, noise, and perturbation is sufficiently 
weak. 

To derive a dynamical equation for the single-oscillator phase PDF of our system, we perform the first, microscopic 
phase reduction of Eq. ([2]). Using the microscopic phase sensitivity of each oscillator, Z((/)) = Vx'/'|x=Xo(0) M' ^^^ 
dynamics of the system can be reduced to the following Langevin phase equation: 



1 ^ 
^jit) = ^ + ]^ E^('^j) • '*^(^o(<^fe)) + V^Z{c^,) ■ r,,(f) +C(0j,t), 



(3) 



fc=i 



where the effect of the perturbation is given by S.{(l)j,t) = Z{(j)j) ■ s(t). From this equation, after averaging and 
taking the limit as A^ ^ 00 (see [2|, |lj| for details), we obtain the nonlinear Fokker-Planck equation (FPE) for the 
single-oscillator phase PDF P{(j),t), 



d d 



w - 



V{^~<j,')P{cj>',t)d<p' + S,{<i^,t) 



P{q^,t)}+D4-:rP{cf>,t), 



(4) 



where the phase coupling function is given by Ticf) — (p') — {1/2tt) J_ Z{(t) + ifj) ■ V{(f)' + ip)dip, and the diffusion 

coefficient by Z? = {Do/27t) J_ Z{ip) ■ Z{ip)dil;. The external perturbation £,{(/>, t) in the drift term is not averaged, 
but kept as a time-dependent term. 

To realize stable macroscopic oscillations, we assume that the coupling function T{(J3) is unimodal and satisfies 
the in-phase (attractive) condition dT{4>)/d(j}\^=o < [4I. When there is no external perturbation (^(0, t) = 0), the 
solution of Eq. ([4]) is given by a rotating wave packet of the form Eq. ([1]), which behaves as follows: In the absence 
of independent noise {Dq — 0), the oscillators are completely synchronized, so that P{(j),t) is delta-peaked. As Dq is 
increased, P{(t>, t) takes on a finite width while rotating at a constant speed. As Dq increases past a critical value Dc, 
the collective oscillation disappears, and P{4',t) = l/27r. For the following analysis, we take < Do < Dc- 

We now derive the collective phase sensitivity C(^) from the microscopic phase sensitivity Z{(t)) by seeking a closed 
equation for the collective phase <&(t). This is performed by the second, macroscopic phase reduction, similarly to 
the derivation of wavefront dynamics for oscillatory media (^. Our treatment here is a generalization of [141] for 
nonlocally-coupled noisy oscillators. Let us ignore the perturbation for the moment {^{(t>,t) = 0). Introducing a 
corotating phase with the wave packet, = (j> — <I>(t) = — fit — O, the solution of Eq. ^, P{(j>,t) = f{0), satisfies 

D{d^/d9^)f{e) - {d/dd) I {uj-Vl) + J^^ T{9 - 9')f{e')d9' f{9)\ = 0. Small deviations u{9) from this unperturbed 
wave packet f{9) obey a linearized equation dtu{9,t) ~ Lu{9,t), where the linear operator L is given by 



M^)-^^"W-^ 



iuj-n)+ / T{9 - 9')f{9')d9' 



u{9) 



±^f{9)fj{9-9')u{9')d9' 



(5) 



Defining the inner product as [u*{9),u{9)] = J_^u*{9)u{9)d9, we introduce an adjoint operator L* of L by 
[u* (9) , Lu{9)] = [L*u*(9),u{9)]. In the calculation below, we only need zero eigenfunctions uo{9) of L and ^0(6*) 



of I/*, which we normalize as [uq{9),Uq{9)] = 1. It is easy to check that uo{9) can be chosen as uo{d) — df{9)/d9, 
reflecting the translational symmetry of Eq. (|4]) [2, [iJI . 

We now incorporate the effect of perturbation, ^{(j),t). To obtain a closed equation for $(i), we assume that the 
solution of Eq. ^ is still given in the form of Eq. ([1]), ignoring its slight deformation at the lowest order. Instead, 
we allow the constant 9 in Eq. ([T|) to vary slowly with time as 6(i). Plugging P{(j),t) — f{9{t)) = fi^cj) — ^t — 9(t)) 
into Eq. g]), we obtain ua{9)Q{t) = {d/d9) {^{9 + 0.1 + Q{t),t)f{9)}. Taking the inner product with u^{9) on both 
sides yields Bit) = J^^ u^{e){d/d9) {^{9 + nt + e{t),t)f{9)} d9. Noting that ^{9 + 01 + e(t), t) = ^{9 + '^{t),t) = 
Z{9 + $(i)) • s{t), we arrive at the macroscopic phase equation obeyed by ^(t), 

m = ^ + c{Ht))-s{t), (6) 

where the collective phase sensitivity is given by 

CW = j\ (-^/(^)) Z{9 + <^)d9. (7) 

Thus, C{^) is expressed as a convolution of the microscopic phase sensitivity Z{9) with a kernel —ul^{9)'f{9). Note 
the kernel is not merely the PDF f{9) as we would naively expect. 

Now let us illustrate our results with an example of N noisy globally-coupled Stuart-Landau (SL) oscillators 
described by the following Langevin equation 0,111]: Wj{t) = (1 + iLOo)Wj ~ {I + iP)\Wj\^Wj + {K/N)Y,k=i Wk + 
y/2Dorij{t) + s{t). Here, the complex amplitude Wj = Wj^ + iWj (R and I express real and imaginary components, 
respectively) describes the state of the jth oscillator (i.e. Xj = {Wj^,Wj)), wq and (3 are osciUator parameters, 
K the coupling strength, s{t) = s^{t) + is^{t) the macroscopic complex perturbation, Dq the noise intensity, and 
rij{t) = ?/?(t) -l- irfjit) the independent complex Gaussian white noise of zero mean and unit intensity which satisfies 
{vf{t)Vkit')) = {v]{t)r]l{t')) = dj,kS{t - t') and {vfit)viit')) = 0- The phase (/> is defined on the complex plane as 
(j) — argW — /3\n\W\, which grows constantly with a natural frequency lu — luo ~ f3 in the absence of coupling and 
other external effects. The microscopic phase sensitivity has both real and imaginary parts, given analytically as 
Z{(f>) — {Z^{(f>), Z^{(j))) = (— sincj) — Pcoscj), coscj) — /3sin0), and the phase coupling function is given by T{(f> — (j)') = 
-K^/T+^sin{(j) - (j)' + a),a = arg(l + if}) 0, [3]. We use the complex order parameter A{t) = i?(t)e**(*) = 
(1/^) Sfe=i s"^'''*'' — II e*'^/((/) — ^{t))d(f> to quantify the collective oscillation, whose modulus R{t) serves as a 
measure of the coherence [3, UM- By choosing /(0) so that J^ e^'^ f{4>)d4> is positive and real, and taking the initial 

collective phase as = 0, namely, P(0, 0) = /(</)), the phase $(i) of the order parameter A{t) matches the collective 
phase $(t) defined in Eq. ([1]). 

Figure [1] shows the stationary PDF /((^), the zero eigenf unctions uo(0) and Uq{(J)), and the kernel —UQ{(j))'f[(j)), 
obtained by numerically solving Eq. ([4]). The parameters ujq = 2, j3 = 1, K = 0.05 are fixed, while the noise intensity 
is varied as Dq = 0.009, 0.011, 0.012. As Dq approaches the critical value Dc = 0.0125 [2, uM, f{,4') becomes 



flat, while the amplitudes of Uq{(I)) and the kernel —UQ{(j))'f{(j)) increase. Figure [21^ a) shows the real part of the 
collective phase sensitivity C^(*5') calculated using the results shown in Fig. [T] (the imaginary part is simply given 
by C{^) = C^(^ — ""/S) due to symmetry of the SL oscillator). For comparison, we show the microscopic phase 
sensitivity, Z^{<l)) with = <I>, which corresponds to the Dq ~> limit. For Do > 0, Cf,($) is different from Z^{(f)) 
due to distributed individual phases. As Dq —>■ Dc, C^(*i?) diverges as (Dc — Do)~^^^ [17|. Figure [UJb) compares the 
theoretical C^(^) at Dq — 0.009 with those directly measured by adding sufficiently weak impulses (|s| — 0.001) to 
the nonlinear FPE ^ and to the Langevin equation Q with N = 10^000 oscillators jl8|. The results of the nonlinear 
FPE corresponding to the limit N ^ oo agrees well with the theory [19|] . The results of the Langevin simulation show 
a wide distribution of the values because N is necessarily finite, though the piecewise average value shows reasonable 
agreement with the theory [20|. 

As an application of the collective phase sensitivity, we analyze the stability of the synchronized state of two 
uncoupled populations of collectively oscillating, globally-coupled SL oscillators due to a macroscopic common noise. 
Consider two uncoupled macroscopic phase oscillators driven by a common, weak Gaussian noise s(i): $('^)(i) = 
0+(^{^^'^\t))-s{t) (cr = 1,2), where the noise correlation is given by (s(i)s(O)) —IC{t). In [l5| , it is shown that a weak 
common Gaussian noise will cause synchronization of uncoupled oscillators, and the Lyapunov exponent that quantifies 
the growth rate of an infinitesimal phase difference is given by A = (l/27r) L dt C{t) J_^ d$ C"(*&) ' C(*& ~ ^t)j where 
A < always holds. For the globally-coupled SL oscillators, the perturbation s{t) has real and imaginary components. 
We use the Ornstein-Uhlenbeck process [16(, z{t) = —z/t + x{t)/''' with xit) a zero-mean Gaussian white noise of 
unit intensity, to create colored Gaussian noises z^'^(t) with the correlation time r, and put s{t) = V2S{z^{t), z^{t)) 
where S controls the intensity of common noise. The correlation function is given as C(i) = 5'exp(— |i|/T)/T. 



Figure [3] shows the results for two uncoupled populations of N — 1000 globally-coupled SL oscillators receiving 
macroscopic common Gaussian noise of intensity S = 0.001 with correlation times r = 1.0, r = 0.5, and t —^ 
(white limit). Other parameters are luq — 2, (3 = 1, K = 0.05, and Dq = 0.007. Figures [3fa) and[3I^d) show the time 
evolution of the real part of the order parameter A(t) of the two populations obtained by the Langevin simulation. At 
time t = 0, the two are not synchronized, but by i = 600, they show near-complete synchronization, demonstrating 
that macroscopic common noise can cause the collective oscillations of two uncoupled populations to synchronize. 
We stress that this is purely a macroscopic phenomenon of the collective phase becoming entrained; the oscillators 
in the two populations never become individually synchronized. Figures Eljb) and ^e) show the distributions of the 
states {Wj} of the microscopic oscillators on a complex plane a,t t — and t = 600, and Figs. [3jc) and[3jf) show the 
corresponding histogram of the phases {(j)j}. At i = 0, the distributions of oscillators only slightly overlap, while by 
t = 600, they overlap considerably. Figure [3Jg) compares the Lyapunov exponent A as measured from the Langevin 
simulation, averaged over 1500 sample paths, with the theoretical results. Despite the relatively large fluctuations 
shown by the collective oscillations, we see that A shows good agreement with theory. As Dq approaches the critical 
value Dc, the amplitude of the collective phase sensitivity C(^) increases, so A takes on ever larger, negative values. 
However, near the critical point, the phase response becomes strongly nonlinear, so that the numerical results diverge 
from the linear theory based on the collective phase sensitivity. 

In summary, we derived the macroscopic collective phase sensitivity from the microscopic phase sensitivity for 
a population of globally-coupled oscillators, and analyzed the common-noise induced synchronization of two such 
populations. By virtue of the assumption that the constituent oscillators are only weakly perturbed, we could utilize 
the phase reduction method to construct a general framework for the collective phase sensitivity, which may bring a 
new perspective to the existing body of research on coupled collective oscillations 10|, |ll|, |l3] . Since we treated the 



SL oscillators as an example, the resulting collective phase sensitivity was always sinusoidal, with only a change in 
amplitude. For other types of oscillators, the shape of the collective phase sensitivity may differ significantly than 
that of its constituent oscillators. Furthermore, the collective phase response to strong macroscopic perturbations 
should prove to be even more intriguing, though our present framework based on phase reduction is not applicable in 
this scenario. More detailed and generalized analysis will be reported in the near future. 
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This is because near a supercritical Hopf bifurcation, the amplitude of the deviation of /((/>) from the uniform PDF l/27r 

scales as (Dc — Do)^ , so that the amplitudes of uq{4>) and the kernel scale accordingly as (Dc — Dq)"^''^. 

An external impulse s changes the phase of each oscillator <j) ^ cfi + Z{<j)) ■ s, so the phase PDF of changes P((jf>) —> 

Z^ P{^)m - (^ + Z(r^) ■ s))d^. 

Because the impulse amplitude |s| is finite though small, there is a slight deviation in the result which becomes more 

pronounced as Do approaches Dc. 

For the discussion of noise-induced synchronization, the mean, and not the fluctuation, is the vital factor. 
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FIG. 1: (color online). Globally-coupled SL model, (a) Stationary PDF /(<;/>), (b) zero eigenfunction Uq{(J)) of L* , (c) zero 
eigenfunction uo{(j)) of L, and (d) kernel —U(){(l))'f {(/)). 
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FIG. 2: (color online). Collective phase sensitivity of globally-coupled SL model, (a) Theoretical curves of (^^(cl?). Single- 
oscillator phase sensitivity Z^'{(j>) is also shown with = $. (b) ^''^($), found from theory and from numerical simulations of 
nonlinear FPE and Langevin equation (LE). For the LE simulation, individual raw responses are shown as dots and piecewise 
mean values (25 bins) are shown as a solid line. 
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FIG. 3; (color online). Common-noise induced synchronization of globally-coupled SL populations. Time evolution of the 
real part of A{i) from (a) i = and (d) f = 600. Snapshots of {W^} at (b) t = and (e) t — 600 (1/5 of population shown). 
PDF of phase at (c) t = and (f) t = 600. (g) Comparison between Lyapunov exponent A, found from theory and Langevin 
simulation. 



